Drag force on a sphere moving towards an anisotropic super-hydrophobic plane 
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We analyze theoretically a high-speed drainage of liquid films squeezed between a hydrophilic 
sphere and a textured super-hydrophobic plane, that contains trapped gas bubbles. A super- 
hydrophobic wall is characterized by parameters L (texture characteristic length), fei and &2 (local 
slip lengths at solid and gas areas), and 4>i and (j)2 (fractions of solid and gas areas). Hydrody- 
namic properties of the plane are fully expressed in terms of the effective slip-length tensor with 
eigenvalues that depend on texture parameters and H (local separation). The effect of effective 
slip is predicted to decrease the force as compared with expected for two hydrophilic surfaces and 
described by the Taylor equation. The presence of additional length scales, L, bi and 62, implies 
that a film drainage can be much richer than in case of a sphere moving towards a hydrophilic plane. 
For a large (compared to L) gap the reduction of the force is small, and for all textures the force 
is similar to expected when a sphere is moving towards a smooth hydrophilic plane that is shifted 
down from the super-hydrophobic wall. The value of this shift is equal to the average of the eigen- 
values of the slip-length tensor. By analyzing striped super-hydrophobic surfaces, we then compute 
the correction to the Taylor equation for an arbitrary gap. We show that at thinner gap the force 
reduction becomes more pronounced, and that it depends strongly on the fraction of the gas area 
and local slip lengths. For small separations we derive an exact equation, which relates a correction 
for effective slip to texture parameters. Our analysis provides a framework for interpreting recent 
force measurements in the presence of super-hydrophobic surface. 

PACS numbers: 83.50.Rp, 68.08.-p, 83.60.Yz 



I. INTRODUCTION 



Super-hydrophobic Cassic (SH) surfaces have opened a 
whole new field of investigation, with both fundamental 
and practical perspectives [l| . They are able to trap air at 
the liquid-solid interface, leading to remarkable ('super') 
properties, such as a very large water contact angle and 
low hysteresis. This strong hydrophobicity has macro- 
scopic implications in the context of self-cleaning [2| and 
impact processes SH surfaces could also revolution- 

ize microfluidic lab-on-a-chip systems [11,0] since the large 
effective slip of SH surfaces |7h3 compared to smooth hy- 
drophobic channels [lol - [T2| can greatly lower the viscous 
dra g. S H surfaces can also amplify electrokinetic pump- 
ing and mixing [^[l^ in microfluidic devices. 

In addition, this superlubrication potential should dra- 
matically modify a squeeze film drainage between sur- 
faces. In the previous study [l7| we have analyzed a drag 
force, F, on a hydrophilic disk approaching a SH wall 
(the Reynolds problem). Here we explore what happens 
when a hydrophilic sphere of radius R is driven towards 
a SH plane with a velocity U (see Fig.[T]), i.e. we address 
the so-called Taylor problem. Beside its significance as a 
geometry of surface forces apparatus (SFA) and atomic 
force microscope (AFM) force experiments, it represents 
a typical situation of phenomena of 'viscous adhesion', 
coagulation, and more. In case of hydrophilic surfaces a 



hydrodynamic force reads |18j 
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when the gap h R. Here ^ denotes a fluid dynamic 
viscosity. In case of a hydrophobic plane, characterized 
by a constant slip length b (the distance between the solid 
at which the flow profile extrapolate to zero) the Taylor 
expression has to be corrected for slip [l^ [lo] 
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The factor /* associated with hydrodynamic slip can 
significantly decrease the hydrodynamic resistance force 
provided h is of the order of 4& or smaller. Eq.(l2]) is of- 
ten used to infer the value of a slip length from the force 
experiment. Indeed, the dynamic SFA/AFM force mea- 
surements [2l[ are extremely accurate at the nanoscale 
as compared to direct flow profiling, or vclocimetry p^ . 
Therefore, it is possible to measure even nanometric slip 
lengths [nlEillil. The advantage of this (hydropho- 
bic vs. hydrophilic) geometry of configuration is that it 
allows to avoid the gas bridging and long-range attrac- 
tive capillary forces [25| , which appear when we deal with 
interactions of two hydrophobic solids [26l - l28| . 
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successfully applied to solve the Taylor problem in the 
presence of the SH surface. Our theory has the merit of 
yielding useful (approximate) analytical results as well as 
being very well suited to numerical work. 

Our paper is arranged as follows: In Sec. II some gen- 
eral consideration concerning a description of a drainage 
of a liquid film confined between a hydrophilic sphere 
and arbitrary SH textured plane are presented. Here we 
also describe some universal asymptotic solutions, valid 
for any texture. Sec. HI contains analytical and numeri- 
cal results for striped SH surface. We conclude in Sec. IV 
with a discussion of our results and their possible rel- 
evance for force experiments. Appendix A contains a 
derivation of equations in a thin gap limit. 



FIG. 1: (Color online) Sketch of a hydrophilic sphere ap- 
proaching a striped super-hydrophobic surface (top), and 
its model representation as a flat interface with patterns of 
boundary conditions (bottom). 



For these reasons it is attractive to consider the hy- 
drodynamic interaction of a hydrophilic sphere with a 
SH surface. However, despite of its importance for force 
experiments and numerous applications, the quantitative 
understanding of the problem is still challenging. The 
heterogeneous nature of the SH texture makes difficult 
a precise discussion of the liquid flow past composite re- 
gions, especially when surface patterns are anisotropic. It 
would seem therefore appropriate to bring a more mod- 
ern theoretical technique to bear on this problem. In this 
paper we present some results of a study of a force, act- 
ing on a sphere approaching to the SH wall. Our theory 
is based on the effective slip approach introduced origi- 
nally to describe flows past a single interface or in a flat 
channel [2^ [s^l . In this approach, the effective slip 
is evaluated by averaging of a flow over the length scale 
of the experimental configuration being applied at the 
hypotetical smooth surface. Such a boundary condition 
mimics the actual one along the true heterogeneous SH 
surface, where gas pockets are stabilized with a rough 
wall texture. A corollary of this is that the effective hy- 
drodynamic slip, boff, gene rally depends on the flow di- 
rection being a tensor [29| and that it also depends on 
the separation between surfaces 0, [T^ . The concept of 
an effective tensorial slip has provided a great deal of in- 
sides into various factors that determines flow in a flat 
channel, and allowed one to obtain very simple solutions 
of complex problems [2^ . We shell see that it can be 



II. GENERAL THEORY 

In this section we describe the theory of a hydrody- 
namic interaction of a sphere with an idealized SH sur- 
face in the Cassic state (sketched in Fig.[T]), where a liquid 
slab lies on top of the surface roughness. The liquid/gas 
interface is assumed to be flat with no meniscus curva- 
ture, so that the modeled super-hydrophobic surface ap- 
pears as a perfectly smooth with a pattern of boundary 
conditions. The latter are taken as low partial slip (6i) 
over solid/liquid areas and as large partial slip (62) over 
gas/liquid regions. We denote as i5 a typical length scale 
of gas/liquid areas. The fraction of solid/liquid areas 
will be denoted 0i = (L — 5)/L, and of gas/liquid area 
(j)2 = 1 ~ 4>i ^ 5 / L. In this idealization, some assump- 
tions may have a possible influence on the friction prop- 
erties and, therefore, a hydrodynamic force. First, by 
assuming flat interface, we have neglected an additional 
mechanism for a dissipation connected with the meniscus 
curvature [sTI - fssj . Second, we ignore a possible transition 
towards impaled (Wenzel) state that can be provoked by 
additional pressure in the liquid phase [13, [l^ . 

The flow of liquid in the gap satisfies Stokes equations 



dz 
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where x-r = Vx^x + VyGy is the lateral velocity, is the 
normal velocity (directed towards the sphere), p(x,y) is 
the local pressure, and Vr is the differential operator in 
a plane (x, y), given by 
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The SH plane (generally anisotropic) exhibits uniform 
tensorial slip [2^, so that the boundary conditions are 



Z = 0: iVr), = (6eff) 



dz 



Vz = 0; (6) 
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z = H{x,y) : v = -Ue^, 



(7) 



with summation over repeated indices. Near the axis the 
surface of a sphere can be approximated by a paraboloid 
of revolution 



H = h + — + Oir^) (r«i?), 
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where = + y^. The solution of Eq.® by imposing 
boundary conditions ^ and ([7]) yields 



(9) 



where i and j ~ {x, y}. 

For convenience, we now diagonalize the effective slip 
length tensor bcff, by aligning x-axis with the 'fast' axis 
of greatest forward slip, which is always perpendic- 
ular to the 'slow' axis of least forward slip, 6^ [l^. By 
integrating Eq.(|3]) with ©-(IT]), and by using we then 
derive a partial differential equation for pressure 



dx 



where fcy _l are the permeabilities of a flat channel of 
thickness H with one SH wall 



Here k = represents an isotropic permeability of 

a flat hydrophilic channel with the same thickness, and 
/cji ^ are corrections to this permeability due to SH slip. 

The solution of Eq. (|T0)) represents a very difficult prob- 
lem since in general broken symmetry of the prob- 
lem caused by a texture anisotropy does not allow one 
standard simplifications, such as representing p{x, y) as 
p{H) dO, H^, etc. Since fc|| j_ are radially symmetric, it 
is more convenient to solve Eq. (|10p by using polar coor- 
dinates 
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with boundary conditions 
p(i?,(^) = 0. 
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The condition for a pressure takes the form p{r ~> oo) = 
0, provided R ^ max{6, L, h}. 

In general case, the solution of Eq.(fTT|) is not radially 
symmetric, Afc ^ 0, due to the terms proportional to 
cos 2lp and sin 2ip. However, it is symmetric with respect 
to X and y axes, hence, it can be presented in terms of 
cosine series: 



P = Pq ('') + Pi cos 2(p + p2 (r) cos Aip + 



(13) 



This expansion can then be used to solve Eq. ([TT]) numer- 
ically. Also an asymptotic solution can be constructed, 
provided e = Ak (0) /(fc) (0) < 1 js^. 

The hydrodynamic resistance force is given by 

oo 27r oo 

F{h) = Ft.!* ^ j j ^"^^^ =2^^^ j Po dH. (14) 

h 

We stress that terms with cos(2n(/?) vanish after integra- 
tion over ip in the interval [0, 27r], so that only isotropic 
part of the pressure, po, contributes to the drag force. 

In two limiting cases, namely small {H <^ L) and large 
{H :$> L) distances between surfaces, the pressure distri- 
bution could be treated with high accuracy as radially 
symmetric [37j . Once such assumption is made, the gov- 
erning equation Ea. (jlip for pressure could be rewritten 
in the form: 



H{k) dp 
r dr 
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Here we use tilde to distinguish approximate p from 
which is an isotropic part of a general solution, Eq.([l^ 
Ea. ([T5)) can be then integrated to give 
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SO that we can represent y) as a function of y) 
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The hydrodynamic resistance force can be then calcu- 
lated as 



F = Ft!* = 2ttR / p dH 



(18) 



In general case, the pressure and the force should be 
found numerically and their calculations require detailed 
knowledge about eigenvalues of effective slip tensor for 
some particular texture. However, in some situations 
asymptotic solutions could be obtained, and they are the 
same for all textures. Below we briefly discuss some of 
them. 
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• For large gap, H ^ L, the corrections to permeabil- 
ity are [H 



kT, ~ 1 
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Even for strongly slipping (62 ^ L) surfaces, the effective 
slip lengths remain independent of the gap thickness and 
are small compared to it [l7j . Consequently, 
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so that Ak/{k) b\/H I. The solution of Eq.dlll) 
can then be constructed in terms of a series of e. In 
this case a non-axisymmetric part of the pressure field is 
small, and the axisymmetric part may be approximated 
by p with a very high accuracy, po = P [l + O (e^)] [13] ■ 
The correction to SH slip is then 
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which coincides with the result, obtained earlier for cor- 
rugated (Wenzel) surfaces via Lorentz reciprocal theo- 
rem [3I]. Physically, Eq. (|20| means that the anisotropic 
SH plane at a distance h from the sphere apex is equiv- 
alent to a no-slip plane at a distance h + (6|)fj + 6^)/2. 
The shift of this equivalent no-slip plane from the real 
SH surface is thus equal to the average of the eigenvalues 
of the effective slip-length tensor. 

For weakly slipping surfaces, 62 ^ L, Eq. ([20l) can be 
further simplified since in this situation the effective slip- 
length tensor is isotropic, and the effective slip coincides 
with the surface average [13, • Therefore, 
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• For thin gap, H ^ L, and weakly slipping surfaces, 
b2 <^ H, the flow is also isotropic [^, and corrections 
to permeabilities are given by Eg. p^ . In what follows 
corrections for SH slip arc given by Eg. (f21]) . In other 
situations /* is sensitive to the pattern geometry and 
should be calculated individually for a particular texture 
of interest. 



hi and 62, vary only in one direction, thus allowing us 
to highlight effects of anisotropy. The problem of flow 
past striped SH surfaces has previously been studied in a 
context of a reduction of pre ssure-driven forward flow in 
thick [39l - l4l| and thin [sl. |42| channels, and it is directly 
relevant for a passive mixing fsl- E^- and a generation of 
a tensorial electro-osmotic flow M, UM ■ 



A. Arbitrary gap, non-slipping solid regions 

The striped SH surface with 61 = currently seems to 
be only one where the eigenvalues of the slip length tensor 
have been calculated for an arbitrary thickness of the 
channel [1, . For the general case the eigenvalues can 
be calculated semi-analytically, but accurate analytical 
results have been found in the thin and thick channel 
limits. 

In the case of thick channels, H ^ L, the eigenvalues 
of the slip-length tensor read [33] 
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Flow in a large channel does not depend on H, and is 
controlled by the ratio of the local slip length b to texture 
period L. When 62/i < 1, Eqs.dm and (gSl) predict the 
area-averaged isotropic slip length, b^^^^ ~ (t>2b2- When 
b/L ^ 1, expressions and take form 
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that coincides with results obtained for the perfect slip 
(62 — 00) stripes [40| . 

In the case of thin channels, H L, striped surfaces 
have been shown to provide rigorous upper and lower 
Wiener bounds on the effective slip over all possible two- 
phase patterns [1] 
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III. STRIPED SUPER-HYDROPHOBIC 
SURFACES 



At b/H S> 1 these give truly tensorial anisotropic effec- 
tive slip 



To illustrate the predictions of the general theory, 
we now focus on a flat, periodic, striped SH surface. 
This canonical texture is convenient to explore the basic 
physics of the system since the local (scalar) slip lengths. 
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but at b/H -C 1 it leads again to a surface average slip. 
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FIG. 2: A correction to the resistance force for SH slip, /*, 
as a function of h/L, calculated with foi = and 62/^ = 10 
(dash-dots), 1 (solid curve), and 0.1 (dashed curve). The 
fraction of gas regions is 02 = 0.5 
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FIG. 3: A correction for SH slip to the resistance force exerted 
on a sphere moving towards the anisotropic plane of stripes 
with no slip on solid- liquid regions (61 = 0). 02 = 0.3 (dash- 
dotted), 0.5 (solid), 0.8 (dash), 62 /-L = 10. 



Calculations of bcs/H, made for several 62/^, were 
used to compute pressure, and then the correction for 
effective slip, /*, as a function of the gap (Fig. [2|). For 
this numerical example we have chosen 02 = 0.5, which 
gives a maximum transverse flow in a thin channel situa- 
tion [1^ . All calculations were performed in the assump- 
tion of the radially symmetric pressure field, p = p{H), 
which should give reasonably accurate results even at 
the intermediate distances, i.e. when the effective slip 
changes rapidly with H [3, H^. Fig. [3] shows that at 
large distances all curves converge to /* = 1, i.e. the 
drag force is the same that it would be in case of a hy- 
drophilic plane, F = Ft- This conclusion directly follows 
from Eq. (PU)) and is valid for any, however large, &2- The 
correction for SH slip significantly decreases when h be- 
comes of the order of L and smaller. It can be seen that 
the increase in 62 leads to a smaller drag force at large 
and intermediate distances. However, at h <^ L all curves 
tends to a constant, which does not depend on the value 
of the local slip length at the gas sectors. 

Fig. [3] includes theoretical curves calculated for 62/^ = 
10 with different fractions of the surface gas phase. Re- 
sults show that the correction for effective slip has a ten- 
dency to decrease with (j)2, and that for each 02 there 
is some minimum value of /* that all curves approach 
asymptotically. We are now on a position to quantify 
this important result. When 61 = and h <^ min{&2ji} 
corrections to permeabilities take a form @ 
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i.e. fcy/A;^ is independent of H, and p coincides with the 
exact solution of Eq. ([TT|) . Therefore, 
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FIG. 4: A correction for SH shp as a function of a fraction of 
the gas phase at &i = and fe/i ~ 10. From top to bottom 
curves with symbols correspond to h/L = 0.01, 0.1, 1, and 10. 
Solid line without symbols shows calculation results obtained 
in the limit of a thin gap, h <^ min{b2, L}. 
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This result coincides with the expression obtained earlier 
for an interaction of a disk with a similar SH surface [l7j 
and shows that in the limit of thin gap /* depends only 
on the surface fraction of the gas phase. 

Since 02 is one of the key parameter determining re- 
duction of drag, in Fig. |4] we plot the correction for SH 
slip as a function of 02 at several h/L. We see that when 
02 is very small, the correction for SH slip tends to its ab- 
solute maximum, /* = 1. In the most interesting limit, 
02 — )■ 1, we can reach the minimum possible value of 
correction for SH slip, /* 1/4, provided 62/^ is large. 
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B. Thin gap, slipping solid regions 

In previous subsection we ignored the (small) slip 
length at the solid areas. The finite bi cannot modify 
our conclusions made for a thick or intermediate gap. It 
could, however, influence a thin gap situation. Fortu- 
nately, in the limit of a thin channel, the exact values 
of longitudinal and transverse permeabilities are known 
for any piecewise constant local slip length, bi and 62- 
These values represent bounds that constrain the attain- 
able effective permeabilities and allows us to calculate 
effective slip lengths Q, and then /* for arbitrary local 
slip lengths. 

The expression for a longitudinal permeability can be 
presented as Q 

(30) 
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Transverse stripes, in turn, satisfy the equation 



where 
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These expressions allow one to calculate effective slip 
lengths in eigendirections 
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In general case, the pressure and the force should 
be calculated numerically with the procedure described 
above. Some special cases, however, allow one exact an- 
alytical solutions. In particular, when H <^ min{62,i} 
the corrections to permeabilities due to SH effective slip 
become 
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The expressions for pressure and force in this situa- 
tion arc rigorously derived in Appendix \^ and an exact 
analytical expression for /* is given by Ea. (jA2[) . Fig. [5] 
includes theoretical curves, /* vs. 02, calculated for sev- 
eral bi/h. At 02 — > 0, and 02 — > 1, our results are consis- 
tent with earlier predictions for homogeneous hydropho- 
bic solids, with bi and &2, described by Eq.© [l^. In 
general, a finite hydrophobic slippage at the solid sectors 
can significantly reduce /* at short distances. 
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FIG. 5: The correction for SH slip, /*, vs. surface fraction of 
tlie gas phase, (j!)2, in the limit of a thin gap and for 62 /fe 3> 1. 
From top to bottom b\/h = 0, 0.5, 1 and 2. 



TABLE I: Some useful asymptotic expressions for the correc- 
tion factor, /*, in case of a striped super- hydrophobic surface 



Case 



h ^ min{62, L} 



&i = 0, <C min{62, L} 



L > /i > 62 



h:$> L 
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with Eqs.jlll) and ([23 



IV. FINAL REMARKS 

Certain aspects of our work warrant further com- 
ments. We have presented data describing the squeeze- 
film drainage of liquid in thin films between a hydrophilic 
sphere and a SH plane. The predicted drag force is 
smaller than expected for two hydrophilic solids and 
shows much more complex behavior as compared with 
expected in case of a uniform hydrophobic wall. For con- 
venience, we now summarize our main results for striped 
SH surface in Table ID At first sight it is somewhat sur- 
prising that a system with one SH surface should exhibit 
such rich hydrodynamic properties. However, it becomes 
almost self-evident (from the presence of extra fundamen- 
tal length scales, such as L, and two local slip lengths) 
that it should exhibit more rich behavior than a geom- 
etry with hydrophilic or even homogeneous hydrophobic 
surfaces. Do our results have any consequence for force 
experiments with real SH surfaces? 

Hydrodynamic force can be measured using various 
techniques, such as AFM, SEA, and more. The geom- 
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etry of configuration is equivalent to a sphere, of a large 
radius of curvature, approaching a flat surface. Hydrody- 
namic forces at small (compared to R) separations have 
been already reported for textured isotropic periodic [U 



and random [2J| SH surfaces. Also we are not aware 



of any measurements of the hydrodynamic force for the 
striped SH surface, measurements with isotropic textures 
lend some support to the picture of the significant reduc- 
tion of a hydrodynamic force presented here. It would 
appear that reduction of a hydrodynamic drag that is 
observed is consistent with an earlier theoretical model, 
Eq.dl]), that gives b^s of the order of 50-200 nm [13, Si, 
which is only slightly larg er than observed at smooth hy- 
drophobic solids [lO, llll I23I I . We suggest that further 
analysis of these measurements or similar measurements 
with other SH textures should employ one of the Eqs. 
given in Table U rather than Eq.(l2]) since the latter makes 
assumptions of a constant, isotropic, homogeneous, and 
independent on distance slip which are not generally valid 
for a heterogeneous SH surface. While we maintain that 
Eq.© is the appropriate formula for describing /* at 
very small (02 0) and very large (02 1) fractions of 
the gas phase we note that other formulas will be required 
in all different situations with the regime of validity given 
in Table n 

From our calculations it is evident that simple far-field 
Eg.ipU)) gives an accurate estimate of /*, and should cer- 
tainly be used to deduce the average of eigenvalues of 
boff at large separations. We note that similar evalua- 
tions could be used for rough or porous surfaces since 
at large distances from the wall the boundary condition 
at the rough interface or fluid- porous interface may be 
approximated by a slip model jH, I45l - l47j . It should be 
possible to make such an analysis of recent data obtained 
with grooved surfaces [48| . which we suspect will signif- 



icantly alter the conclusion of these authors. Note that 
the average eigenvalues of an effective slip length tensor 
at large (compared to L) separations should coincide with 
the measured by velocimetry [i^ or other far-field meth- 
ods We remark and stress, however, that the near- 
field effective slip often measured in AFM / SFA should be 
much smaller than a far-field slip. This is exactly what 
has been observed in experiment [23, [4^ . 
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Appendix A: Analytic solution for p and /* in the 
thin-channel limit 



In this Appendix we derive analytical expressions for 
a dimensionless pressure distribution P = ph'^ {3iJ,UR)~^ 
and for /* in a situation when h min{&2, L} and the 
corrections to effective permeabilities are given by Eqs. 

((Ml), dSzj). 

We use the substitution into Eq. (fT7|): 
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(3 — x) (4 — (j)2x) dx 



0102 (4 ~xf + A (02 + 02 + 1) (4 _ a;) + 160102 



X 



1 



3 , (4-30i) 
2 



2x 



+ (3^ci In 



0102X2 _ _,_ 32 



+ /3^C2 In 



4i02X — 4 — d 



X - 4 + 



1 , -, 1 1]'^ + air] + a2 
Cl in 



C2 In 



V + aa 
77 + 04 



01 (301 - 5) 01 (301 - 1) (201 - 3) 

, C2 = 



32/32 



40/3 32/32 

, 02 = , 03 

05 a5 



8/32a6 
/3 (ae - 4) 



/3(«6+4) 

30102 -06-4' 



30102 -1-06—4' 
05 = 8 + 902 - 90^, 06 = 4a/i - 202 + 202. 



(Al) 



8 



Then one can also integrate analytically Eq. (IT4l) with P given by (jAip to get 



Pdi] 



3 / 1 ^ 1 ^ 



?7 + — ^ ) In [if + aiTj + 02) + a7 In 



2r/ + Qi + 207 
2?7 + ai — 207 



+ C2 [(?7 + as) In (77 + as) - (?/ + 04) In (?/ + 04)]}!^^ 



1 + ^ ) In (1 + ai + 02) - ai + [cia7 + C2 (1 + 04)] In (1 + 04) 



[cia7 + C2 (1 + 03)] In (1 + as) - C2 (04 - as)} , 



— ai — 277 In 77 



(A2) 



07 = 3/3a6/a5. 

I 



Eq. reduces to Eq. ([2]) as 02 ^ and to Eq. ([29ll 
as 0. 
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Supplementary information. 
Drag force on a sphere moving towards an anisotropic super-hydrophobic plane 



Evgcny S. Asniolov, Alckscy V. Bclyacv, and Olga I. Vinogradova 



NUMERICAL AND ASYMPTOTIC SOLUTIONS 
OF EQUATION FOR PRESSURE 



Eq.(ll) of our main paper has the foUowing form: 

H{k) dp H{k) d^p 



dr 



(1) 



+ cos 2ip 
sin2iy9 



or V or 



2HAk 



d'^p 



HAkdp HAkd^p' 
r dr dip'^ 
d{HAk) 2HAk\ dp 
r J dip 



drdtp \ dr 
and its solution ean be presented in terms of cosine series: 



P = Po (f ) + Pi (r) cos 2ip + p2 (r) cos Aip + 



(2) 



Here we derive equations for functions p„ (r) in series 
Eq.([2]) and describe their numerical solution. When Ak 
is small, we construct an asymptotic solution. 

To solve Eq.([T]) we substitute ^ into ([T]) and collect 
terms proportional to cos 2nip. To eliminate singularities, 
we then use logarithmic substitution for variable r and 
introduce the following dimensionless quantities: 



In 



/hR 



Pr. 



TJ^k) 



H 
'h 



D (,7) = - 



exp (20 



77^ Ak 
Afc(O) 



This leads to a system of ordinary differential equations 



Lo-Fb + ^Lq Pi — 



Ak 



(k) (0) 



exp (20 : 



(3) 



L„P„ + eL+Pn+i + eL„ P^-i =0 for n > 0, (4) 

which is expressed in a more compact form by using dif- 
ferential operators 



d fD d\ , ^ 
+ (n + l)^ + 2n {n + l)D, 



d f D d\ ^ d 

-{n-l)^ + 2n {n-l)D, 
etc; 

To solve a finite-difference version of ODE system 
([3]), (jH we resolve the truncated system with respect 



to d^Pn/d^^ 



1,..., iV, by using the Gauss routine, 



and then integrate numerically the obtained system by 
using the fourth-order Runge-Kutta method. This sys- 
tem has a boundary condition P„ — > both at small and 
large r, or at ^ — >■ —00 and ^ — >■ +00. Therefore, the 
dimensionless gap and the permeabilities take a form 



ry ^ 1, dG/d^ = dD/d^ 



as 



-00, 



exp (20 



G, D-jf ~ exp (60 as ^ 



-00. 



Asymptotic linearly independent solutions of ((3)) , (HJ can 
be then foimd in the form: P„i = a^j exp (A^O 7 * = 
1, 27V as ^ ^ -00 and P^ = a^exp(A°°0 as ^ ^ 
+00. The eigenvalues, A° and A°°, and eigenvectors, a° 
and a°°, were found using IMSL routine DEVCRG. The 
system admits both exponentially growing and decaying 
solutions Pni whereas the boundary conditions require 
the decaying one. Thus, we choose solutions with A° > 
and X°° < 0. To resolve all linearly independent solutions 
properly at |^| ^ 1 the orthonormalization method [H, 
is applied. The boundaries of the integration domain are 
set at Tniin ~ 0.01 and J'max = 20, the number of taken 
into account harmonics is = 4. The calculations show 
that the expansion converges fast with A^, since the ratio 
\Pn+i/Pn\ is usually small for all n, which is illustrated 
in Fig. m 

The solution of Eq. ([T]) is simple if ratio k\^/k± (and 
hence Ak/{k)) is constant, i.e., when the permeabilities 
depend on r similarly, fcj| , k*^ ^ g [r) . This is the case in 
a thin channel with bi =^0, when all the permeabilities 
arc proportional to [3]. Then the solution includes 
only axisymmctrie term pq. It can be verified directly 
that solution 



dp 

dr 



IxUr 



12tiUr 



2^(fc> (^k*{H) + kl{H] 



PjH) 
UnUR 



dH' 



H 



'3 {h 



(5) 



(6) 



2 



0.25 




-0.25 



r(hR) 



FIG. 1: Functions P„ (r) = pnh^/S^iUR in series ([2| for 
h ^ min{62, 1/}, 6i/ft = 0.5, 02 = 0.8. Solid line corresponds 
to an analytical solution for P = ph? /'i^U R, circles - to a 
numerical evaluation of Po, dashed line - to 10'^ Pi and dash- 
dotted hue - to 10^ P2. 



to the drag force, differs from the approximate solution 
P given by Eq. ([6]) by the value O {e^) . Substituting 
the expansions into the system ([3]), (jH) and collecting the 
terms proportional to one obtains equations governing 
the first three terms of expansion: 



4fc 



(fc) (0) 



exp (20 : 



=0 as n > 0, 



P^^' =0 as n 7^ 1, 



T p(2) _ T + P'- 



(1) 



(7) 



satisfies Eq. ([T]) since the boundary conditionsare homo- 
geneous, dp/d(p — 0, and terms 



d 

dr \ dr 



Ak 

m 



Mi/, 



HAk dp Ak 



dr 



^lU, 



cancel out. The analytical expressions for P = 
ph^{3^UR)~^ and corresponding resistance force in the 
case when h <C min{62, L} are obtained in the Appendix 
section of the main paper. 

The first three harmonics in series ([2]), plotted as func- 
tions of r/^/hR, are presented in Fig.[TJ It shows that the 
axisymmetric part of pressure distribution is very close to 
the approximate solution P while the non-axisymmetric 
part is several orders less. 

The permeability difference is typically small, |AA:| ^ 
|(fc)|, for example, for large distances between surfaces, 
H ^ L. The asymptotic solution of Eqs. (01) can be 
found in this important case as power series expansion 



Pn 



J2j=o^'' '^ith respect to small parameter e. 
We do not construct the full asymptotic solution, but 
show that the isotropic part Po, which only contributes 



L2P2 



(2) 



'^2 Pi 



(1) 



Pp' = 0, P,[-'> =0 as n>2. 

The equation ([7]) is the dimensionless form of Eq. ([5|). 
For further terms, one can deduce that 

Thus the first-order correction to the pressure field is 
ePi^\ while the correction to the isotropic part e^Pq^' is 
much smaller. 



J(2) _ 
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